Cone beam CT-based adaptive intensity modulated proton therapy assessment using automated planning for head-and-neck cancer

Background To assess the feasibility of CBCT-based adaptive intensity modulated proton therapy (IMPT) using automated planning for treatment of head and neck (HN) cancers. Methods Twenty HN cancer patients who received radiotherapy and had pretreatment CBCTs were included in this study. Initial IMPT plans were created using automated planning software for all patients. Synthetic CTs (sCT) were then created by deforming the planning CT (pCT) to the pretreatment CBCTs. To assess dose calculation accuracy on sCTs, repeat CTs (rCTs) were deformed to the pretreatment CBCT obtained on the same day to create deformed rCT (rCTdef), serving as gold standard. The dose recalculated on sCT and on rCTdef were compared by using Gamma analysis. The accuracy of DIR generated contours was also assessed. To explore the potential benefits of adaptive IMPT, two sets of plans were created for each patient, a non-adapted IMPT plan and an adapted IMPT plan calculated on weekly sCT images. The weekly doses for non-adaptive and adaptive IMPT plans were accumulated on the pCT, and the accumulated dosimetric parameters of two sets were compared. Results Gamma analysis of the dose recalculated on sCT and rCTdef resulted in a passing rate of 97.9% ± 1.7% using 3 mm/3% criteria. With the physician-corrected contours on the sCT, the dose deviation range of using sCT to estimate mean dose for the most organ at risk (OARs) can be reduced to (− 2.37%, 2.19%) as compared to rCTdef, while for V95 of primary or secondary CTVs, the deviation can be controlled within (− 1.09%, 0.29%). Comparison of the accumulated doses from the adaptive planning against the non-adaptive plans reduced mean dose to constrictors (− 1.42 Gy ± 2.79 Gy) and larynx (− 2.58 Gy ± 3.09 Gy). The reductions result in statistically significant reductions in the normal tissue complication probability (NTCP) of larynx edema by 7.52% ± 13.59%. 4.5% of primary CTVs, 4.1% of secondary CTVs, and 26.8% tertiary CTVs didn’t meet the V95 > 95% constraint on non-adapted IMPT plans. All adaptive plans were able to meet the coverage constraint. Conclusion sCTs can be a useful tool for accurate proton dose calculation. Adaptive IMPT resulted in better CTV coverage, OAR sparing and lower NTCP for some OARs as compared with non-adaptive IMPT.


Background
Radiation therapy can serve as a primary or adjuvant treatment to surgery for head-and-neck (HN) cancer patients.Both intensity-modulated radiation therapy (IMRT) and intensity-modulated proton therapy (IMPT) are well suited for the complex anatomy in HN cancer treatments as they deliver highly conformal dose to tumor volumes while sparing organs at risk (OAR).IMRT is a common treatment technique, but studies have demonstrated that IMPT can deliver a superior dose distribution compared to IMRT due to the physical property of proton dose deposition known as the "Bragg peak" [1,2].However, IMPT dose delivery is very sensitive to anatomical variations, setup deviations, and range uncertainties due to Hounsfield unit (HU) to stopping power (SP) conversion.These geometric and positional changes can result in OAR toxicity and under dosage to the clinical target volumes (CTV).Additionally, significant inter-fractional anatomic changes due to weight loss commonly occur during HN treatment [3][4][5].One technique for managing these changes is to adapt the treatment plan to the new anatomy so an optimal dose distribution can be maintained throughout the treatment course.
Traditional offline adaptive radiation therapy (ART) typically requires an acquisition of repeat CT (rCT) scans during the treatment course.Acquisition of rCTs contributes extra imaging dose to the patient, requires additional departmental resources, and might provide a false indication for adaptation if the patient setup is not done carefully [6].In image-guided radiation therapy (IGRT), pretreatment cone-beam computed tomography (CBCT) images are becoming a standard of care for patient setup.CBCT images acquired during the standard care path require fewer resources and reduce the overall imaging dose burden to the patient compared to acquiring additional rCTs.However, it is not recommended to use CBCTs directly for dose calculations due to poor image quality caused by increased scatter, motion, beam hardening, and other imaging artifacts, especially for proton therapy [7][8][9][10].Despite these limitations, early studies exploring the use of CBCTs for accurate proton dose calculations show promise [7][8][9][10][11][12][13][14][15][16][17][18][19].
Attempts have been made to reduce the uncertainties in using CBCTs for ART.Scatter correction on CBCT has been studied by several groups as a reliable method for accurate proton dose calculation [7,9,10].Another approach is using deformable image registration (DIR) to transfer the HU information from the planning CT (pCT) to the CBCT to create a synthetic CT (sCT) [10][11][12][13][14][15].Work by Kurz et al. compared proton doses recalculated on sCT and scatter corrected CBCT and found high dosimetric agreement for both head and neck and prostate treatments [10].
The potential benefits of adaptive proton therapy (APT) have been reported by several groups.Simone et al. compared the APT and non-APT using rCT for HN patients and indicated that the APT improves OAR sparing [1].Another study by Gora et al. compared dosimetric benefits between photon ART and APT for six HN patients and found that APT plans improved target coverage and OAR sparing for brainstem and spinal cord [25].Lalonde et al. compared robustly optimized IMPT plans to daily adaptive IMPT without robustness constraints for ten HN patients and concluded that daily adaptation resulted in better target coverage and OAR sparing [22].Botas et al. investigated the feasibility of a fast weight-tune online adaption approach based on Monte Carlo methods with CBCT and indicated significant improvements in plan quality [26].Nenoff et al. studied the benefits of APT by simulating different nasal cavity filling and setup scenarios for five paranasal patients, concluding that APT improves plan robustness [21].More recently, an investigation by Borderías-Villarroel et al. explored the dose difference between non-APT and APT based on different strategies (manually full re-optimization or automatic isodose volume dose restoration) for ten HN patients.Results demonstrated that the dose restoration method could achieve CTV coverage for half of the patients, but manual full re-optimization was still required for the other patients [23].Bobic et al. compared the daily adaptation and weekly adaptation using another dose restoration method and recommended weekly adaptation achieved satisfactory CTV coverage for most patients [27].While these studies have demonstrated the potential dosimetric benefits of APT with small sample sizes, the potential clinical relevance of these dosimetric benefits has not been explored.
Daily online dose evaluation and plan adaptation requires a fast turnaround to be clinically feasible.In an ART workflow, the contouring and re-optimization steps present a challenge as they are resource intensive, but DIR presents an option to quickly transform contours from the pCT to the daily image.Previous studies have demonstrated the feasibility of using DIR-based sCT for proton dose calculation in HN region via stopping power (SP) comparison, water equivalent thickness (WET) comparison, or gamma analysis between sCT and reference CT [10,[12][13][14][15], but only a few studies investigated the impact of uncertainty of DIR-propagated contours on the dose evaluation or plan adaptation in HN region [15,20].
In terms of the proton plan adaptation, previous studies investigated re-optimization with the same objective list [21,22], manual full re-optimization [1,23], or the fast dose restoration method [23,24] to generate a new plan on the daily image.However, manual full re-optimization is time consuming and dose restoration or re-optimization with initial objective list might not result in the optimal dose distribution when large anatomical changes occur [23].Our previous study investigating an automated planning software model for efficient HN IMPT plan generation demonstrated the model-generated plans have a quality that is, at minimum, comparable to the manual plans produced by dosimetrists [28].
The primary goals of this study are to quantify the accuracy of proton dose calculation using DIR-based sCT and the impact of DIR propagated contour inaccuracy on dosimetric evaluation.We also aim to quantify the dosimetric benefits of an efficient offline adaptive IMPT workflow using weekly CBCT and a pre-validated automated planning software [28].This study utilizes the automated planning system for plan adaptation to ensure the optimal dose distribution and reduce the inter-operator variation during re-optimization.The clinical significance of the dosimetric benefits will be estimated by employing radiobiological modeling.

Patient cohort and IMPT planning
Twenty patients with advanced HN cancer previously treated with IMPT or VMAT were included in this study.All patients were enrolled in a retrospective institutional review board (IRB) approved protocol.For each patient, contrast and non-contrast pCTs were acquired on the same day in a supine position with 2-mm slice thickness using the Siemens Somatom CT 16 slice or 64 slice simulators (Siemens Healthineers AG, Germany).All gross tumor volumes (GTVs), clinical target volumes (CTVs), and OARs (spinal cord, brainstem, parotids, constrictors, mandible, cochlea, larynx, carotids, oral cavity, submandibular, eyes, and optic nerves) were delineated on the contrast CT.The volumes were then rigidly transferred to the non-contrast CT.Patients were prescribed one to three dose levels ranging from 56 to 70 Gy delivered in 30-35 fractions.Table 1 shows the characteristics of the patients included in this study.The CTVs were located in the mid/lower neck area for a vast majority of patients in this study.A daily pretreatment CBCT scan for setup was obtained with a ProBeam compact or TrueBeam onboard CBCT imager (Varian Medical System, Inc, Palo Alto, California).Ten patients had at least one rCT scan using the same protocol as their pCT during their treatment course.[28].The beam number and arrangements were selected based on tumor anatomy and location and field-specific targets were created for each field.Proton spots could only be placed in the field-specific targets that were created to encompass all CTVs including 3 mm positional setup uncertainty and 3% range uncertainty for each field [29].The field-specific targets were modified to prevent beams entering through the chin and teeth area to reduce proton range uncertainty from the movement of the chin or tongue.Streaking artifacts caused by dental implants were delineated and overridden to an HU value approximate to the surrounding soft tissue.The non-linear universal proton optimizer (NUPO 16.02, Eclipse, Varian Medical Systems) was utilized for optimization along with the proton convolution superposition algorithm for dose calculation (PCS 16.02, Eclipse, Varian Medical Systems).A 2 mm × 2 mm × 2 mm dose grid was used along with a relative biological effectiveness (RBE) of 1.1 to weight the dose.All IMPT plans were robustly optimized using ± 3 mm setup uncertainty (in cardinal directions) along with ± 3% proton range uncertainty, resulting in 12 uncertainty scenarios.The targets were the only structures selected to be robustly optimized.Plans were optimized using an autogenerated objective list by the RPP model.If needed, one to two additional optimization iterations were performed to improve the CTV coverage and OAR sparing using a fine-tune of the objective list.These additional optimizations for some patients included satisfying the CTV coverage, reducing target maximum dose, or slightly improving OAR sparing if the clinical constraint was not met.All IMPT plans were normalized such that 95% of the primary CTV volume was covered by 100% of the prescription dose (V100 = 95%).Robust evaluation was performed by introducing the same uncertainty as in the robust optimization.The worst-case scenario in the robust evaluation required all CTVs to achieve at least 95% of each volume receiving 95% of the prescribed dose (V95 > 95%).The dose-volume constraints for OARs are given in Table 2.

Validation of proton dose calculation based on sCT
Thirteen rCTs from ten patients were used to validate the accuracy of proton dose calculation on the sCT.The workflow of the validation is shown in Fig. 1.All the data were imported into an imaging software package (Velocity, version 4.1, Varian Medical Systems, Palo Alto, CA) that contains a B-spline based DIR with a built-in CBCT correction algorithm.For sCT generation, the daily CBCT acquired on the same day as the rCT was rigidly registered to pCT by applying the shifts from the patient setup.A displacement vector field (DVF) was then calculated through DIR between the pCT and daily CBCT.The DVF was used to create a sCT by transferring the HU and overridden HU values (e.g., dental artifact) from the pCT to the CBCT frame of reference.Areas outside the CBCT field-of-view (FOV) were filled with the corresponding image data from the pCT source image.
Visual assessment of the DVF was performed to ensure that the DIR transformations were physically and anatomically reasonable.Anatomical landmarks including bony anatomy and the body contour were confirmed to match in the images.Additional refinement was rarely needed to improve DIR accuracy, and this was accomplished by adjusting the region of interest to focus on a smaller region around the CTVs.Contours from the pCT were transformed to the sCT using the same DVF, as well.Although the DVFs appeared visually accurate, additional corrections to the DIR-propagated contours might be required by a physician based on their review.In these cases, two contour sets were created on the sCT-purely DIR propagated contours and physician-corrected contours.The physician-corrected contours served as the gold standard.Finally, the rCT was deformably registered to the daily pretreatment CBCT.The resulting deformed rCT (rCT def ) served as the gold standard for proton dose calculation as this image provided true HU values in the same frame of reference as the sCT.The gold standard physician-corrected contours from the sCT were rigidly transferred to the rCT def .
The initial IMPT plan was recalculated on both the sCT and rCT def .To validate the accuracy of proton dose calculation on the sCT, 3D global gamma analysis [30] with 3 mm/3% and 2 mm/2% was performed within the region receiving > 10% of the prescribed dose.The gamma analysis was performed in MATLAB R2022b (MathWorks Inc.Natick, MA).

Validation of DIR contour propagation
To investigate the impact of DIR propagated contour inaccuracy (including SP transfer uncertainty) on dosimetric evaluation, the dose from the initial IMPT plan was also calculated on the sCT with uncorrected DIR propagated contours (sCT_uncorr).The dose volume indices from the sCT_uncorr were then compared to the indices from the gold standard rCT def with physiciancorrected contours.Additionally, the dice similarity coefficient (DSC) was calculated for the CTVs and OARs to quantify the agreement between the DIR propagated (uncorrected) and physician-corrected volumes.

Benefit of adaptation versus non-adaptation
Twenty patients with daily CBCTs were included in the assessment of offline APT versus non-APT.Previously published work demonstrated that weekly CBCTs accurately estimate and represent the overall delivered dose to HN cancer patients [3].Therefore, daily pretreatment CBCTs were selected every five fractions (e.g., fraction 1, 6, 11, 16, etc.) to serve as weekly images and were used to create weekly sCTs.
For the non-adaptation (non-adapt) group, the proton doses were recalculated on the weekly sCTs using the initial IMPT plan.The dose accumulation was performed by warping the weekly recalculated dose to the pCT through the inverse DVF that was used to create the sCTs.For the adapted (Adapt) group, the IMPT plans were reoptimized on each weekly sCT using the automated planning software RPP [28].Plans were optimized with the same robustness parameters used in the initial plan and dose accumulation was performed using the same method as for the non-adapt group.The weekly dosimetric parameters for CTV and OARs were compared between the non-adapt and adapt groups.The percentage of fractions that met the dose constraints from Table 2 were calculated for each group.Next, the accumulated dose-volume indices from non-Adapt and Adapt patients were evaluated.Finally, normal tissue complication probability (NTCP) models were employed to estimate the probability of > grade 2 larynx edema, > grade 2 dysphagia, > grade 4 xerostomia, and > grade 2 acute esophagitis following the end of the treatment course [31][32][33][34].All statistical analysis was performed using a two-sided paired t-test in JMP

Validation of proton dose evaluation based on sCT
An example case with different pCT, rCT, CBCT, rCT def and sCT is provided in Fig. 2. High dosimetric agreement was observed between the dose calculated on the rCT def (gold standard) and sCT with a mean gamma pass rate of 97.9% ± 1.7% for 3 mm/3% criteria and 93.7% ± 4.2% for 2 mm/2% criteria throughout the whole body.
The physician-corrected contours on the sCT reduced the deviation from the gold standard, as expected.Table 3 shows that there was no statistical difference between the gold standard rCT def and the sCT with physician-corrected contours (except for the tertiary CTV as described below).The deviation range of the relative Dmean was (− 0.9%, 2.18%) for the larynx, (− 2.37%, 2.19%) for the constrictor, (− 0.12%, 1.82%) for the Parotid_ips, (− 1.51%, 0.81%) for the Submand_ips, and (− 3.28%, 7.53%) for the ipsilateral cochlea.The maximum deviation of V100 and V95 for both CTV_primary and CTV_ secondary was reduced to within 3% and 1% by using corrected contours.It should be noted that for some patients, part of the CTV_tertiary was out of the CBCT FOV when if it extended inferiorly to shoulder region.In these cases, the images out of the CBCT FOV were directly copied from the pCT to sCT.The deviation range for CTV_tertiary V100 and V95 was (− 6.2%, 0.33%) and (− 3.97%, − 0.02%), even with the corrected contour.As a result, it should be noted that it can be error prone to use the sCT to estimate the coverage of CTV_tertiary when it extends inferiorly.For the maximum dose to each OARs, the sCT without contour correction can produce large deviations for brainstem and optic nerves (ON).After the correction of contours, the deviation was smaller as compared to non-corrected contours for most OAR Dmax except ipsilateral eye.

Adaptation versus non-adaptation
The fluctuation of weekly coverage for the CTVs and OAR dose-volume indices with or without adaptation can be found in Figs. 5 and 6, respectively.Without plan adaptation, 6 of 134 fractions (4.5%) for CTV_primary, 4 of 97 fractions (4.1%) for CTV_secondary and 26 of 97 fractions (26.8%) for CTV_tertiary failed to meet the dose constraints of V95 > 95%.In the adaptive cohort, the coverage constraints were met for all the fractions.The percentage of fractions failing to meet the dose constraint was reduced from 23.3% (30 of 129) to 17.8% (23 of 129) for constrictors and from 19.3% (16 of 83) to 1.2% (1 of 83) for larynx when plan adaptation was applied.
Figure 7 shows the difference in NTCP between adapt and non-adapt groups.It was notable that the NTCP of larynx edema was significantly reduced by 7.52% ± 13.59% in the adapt cohort as compared to the non-adapt cohort, while achieving a maximum reduction of 45% for patient 8 with significant weight loss, as shown in Fig. 7.For other structures, there was no statistically significant difference in NTCP between the adapt and non-adapt groups.
The number of structures included in the analysis are indicated in the first column as n.It should be noted that some structures were excluded from analysis as they were far from CTVs and received a cumulative dose < 1 Gy.The maximum dose (Dmax) and mean dose (Dmean) are reported as relative dose.A P value < 0.05 is highlighted in bold and indicates statistical significance.A P value < 0.1 is highlighted in italic  for constrictor, larynx, and bilateral parotid.A reference line is placed at constraint level for each OAR.Note that the figure does not present the dose-volume indices for some patients when the dose was too low (< 1 Gy) or when the parotid or larynx was the primary target An example case (Patient 3) with a large tumor growth is illustrated in Fig. 8.With the large tumor growth, the non-Adapt accumulated dose showed considerable under dosage to CTV and OARs as compared to the planned dose (Fig. 8a).With the plan adaptation, accumulated dose was comparable to the planned dose, as shown in Fig. 8b.Another example (Patient 8) is shown in Fig. 8c.This patient encountered significant weight loss during treatment course.The Non-Adapt accumulated dose greatly increased the dose to the larynx and constrictors, as shown in Fig. 8c, while the Adapt accumulated dose resulted in reduction in dose to the OARs with a comparable dose distribution to the planning dose (Fig. 8d).The mean dose to the larynx, constrictor, ipsilateral parotid, and contralateral parotid increased by 14.13 Gy, 7.75 Gy, 3.69 Gy and 2.52 Gy, while the plan adaptation reduced the mean dose by 11.08 Gy, 7.31 Gy, 1.02 Gy and 1.35 Gy, respectively.

Discussion
The use of sCTs for proton dose calculation based on various methods has been explored by several groups [10][11][12][13][14][15].Thummerer et al. compared proton dose calculation with corrected CBCT using several methods and concluded that accurate proton dose calculation can be achieved with both DIR-based and deep learning-based sCT in HN region [19].The high gamma passing rate for this study between dose calculated on sCT and the gold standard demonstrates that the DIR-based sCT provides accurate SP information for proton dose calculation, which was consistent with other studies [10][11][12][13][14][15].Although we assumed CBCT and rCT had the same anatomical information, there were some discrepancies in cavity filling or emptying observed between the CBCT and rCT def for some patients due to changes in the tongue position.To mitigate its impact, the field-specific targets prevented beams entering through the chin and teeth area.Lack of a robust automatic contour propagation method was still an obstacle for online adaptation.We found that the DIR-propagated contours achieved acceptable quality as compared to gold standard contours (DCS > 0.8) for most structures while more attention should be paid to constrictors and small structures (e.g., cochlea).The dosimetric difference between rCT def and pCT represented the deviation caused by anatomical changes and setup error, while the difference between rCT def and sCT with corrected contours was due to the SP transfer uncertainty arising from DIR inaccuracy.We found that the dose deviation caused by anatomical changes and setup error was much higher than the deviation caused by DIR inaccuracy.The errors introduced by using sCT to estimate the dose-volume indices can be further reduced by using physician-corrected contours for most structures, which is consistent with the previous findings [15].It was observed that for structures proximal to CTV (e.g., ipsilateral parotid), the contour correction by physician was more impactful to the DVH evaluation.Estimation of the dose-volume indices for small structures (e.g., cochlea), even with corrected contours, was still unreliable.Although using the uncorrected contours on sCT to estimate CTV coverage can potentially yield acceptable accuracy, it was still beneficial to perform contour corrections for the CTVs.The CTV was one of the most crucial structures and required careful attention.Of note, the CTV_tertiary was typically out of CBCT FOV when it extended inferiorly to the shoulder region.In these cases, the corresponding information from the pCT was filled in and could lead to possible overestimations of the CTV coverage.An extended FOV CBCT scan to include the shoulder region is recommended to improve the accuracy of dose evaluation for CTV tertiary based on sCT.It should be acknowledged that relying solely on a simple DVH comparison may overlook significant dose differences, particularly in small volume regions near the end of the treatment range.
A study by Kurz et al. investigated sCT with DIR generated contours for HN replanning and found improvement of plan quality with significant hotspot reduction and partial improvements in OAR sparing [20].Another group investigated the use of uncorrected contours on sCT for replanning of lung patients and indicated that this approach can still restore the CTV coverage and reduce the hotspots [35].Therefore, in the context of using CT images for dose evaluation or replanning, whether with or without corrected contours, it is crucial to have a clear understanding of the associated uncertainties, as highlighted in the study.One of the limitations of this study was the small sample we included for the validation of sCT for proton dose calculation (n = 13).
In the present work we found that the constrictors and larynx significantly benefited from APT.The improved target coverage was a common finding in previous studies looking into the benefits of APT [20-22, 25, 26].Gora et al. explored the benefits using six HN patients and found that the APT reduced the hotspots in brainstem by 8 Gy and spinal cord by 14 Gy [25].The current study included a larger cohort of patients with CTVs in the mid/lower neck area for the majority of patients.However, there were no significant differences between non-Adapt accumulated Dmax and Adapt accumulated Dmax or planned Dmax for spinal cord and brainstem.It was worth noting that the benefits of adaptive proton therapy can be largely patient dependent.In this study, it was observed that patients who underwent weight loss can suffer from extreme overdosage to OARs (e.g., patient 8).
The NTCP models employed reveal that plan adaptation would significantly reduce the probability of larynx edema.30% of the patients would benefit from APT with reduction of larynx edema by more than 5%, including the maximum reduction of 44.97%.For dysphagia, 10% of the patients achieved reduction of NTCP by more than 5% using adaptive proton therapy, with the maximum reduction of 13.33%.For xerostomia and acute esophagitis, adaptation did not provide an NTCP benefit of more than 3% compared to non-Adapt for all patients.For some patients, in order to restore target coverage, plan adaptation increased the OAR dose (e.g., patient 2), and hence an increase in NTCP was acceptable.It should be noted that the NTCP value could be improved if the optimizer was given more flexibility during re-optimization, this can be minimized by employing the automated planning model for objective list generation.
In this study, we applied the setup uncertainty during re-optimization.It might be reasonable to assume there was no setup error for HN patient with online adaption, which would provide further OAR sparing [22].While our previous work demonstrated that the weekly sampling was sufficient to represent the total delivered dose, it should be acknowledged that it was one of the limitations of this study [3].It might be was possible that weekly sampling could miss weight loss occurring in a short period of time.Additionally, dose accumulation error could be feathered out if more fractions were employed [36].
In this work, a validated proton-specific automated planning system was employed for both initial planning and re-optimization.It was shown in our previous study that the trained automated planning model can generate IMPT plans comparable to expert plans [28].Manual full re-optimization might be the best option to regain the optimal dose distribution, but it was only applied for offline adaptation.Another option for adaption was to use the initial objective list for re-optimization, but there was a risk of not achieving the optimal dose distribution with significant anatomical differences.Fast dose restoration methods have been investigated by several groups for online adaptive proton therapy because of its high speed for plan adaptation and has been demonstrated to improve the plan quality compared to the nonadapted treatment [24,26,37].Borderías-Villarroel et al. compared dose restoration and manual full plan optimization for plan adaptation and found that dose restoration avoided full optimization for 52% of patients but the remaining patients with a large anatomical change or inaccurate positioning still needed full offline optimization [23].In this study, we employed an automated planning software for re-optimization first to diminish inter-operator variation during re-optimization and second to generate the optimal dose distribution.Since this was still a full re-optimization, the time required was around 15-20 min.However, compared to the manual re-optimization, which requires iteratively tuning the objective list and demands hours, utilizing the model substantially reduces the workload.sCT generation and contour propagation add an additional 5-10 min.With manual contour adjustment, 5-10 additional minutes were required, depending on the DIR quality and number of critical organs required for manual adjustment.In its current form, this method is feasible for offline adaptation, but application to online adaptive proton therapy still requires additional validation including investigation of a fast online QA procedure.Nevertheless, we did observe the benefits of employing this model for restoration of CTV coverage and OAR sparing.

Conclusions
This study assessed the feasibility and potential benefits of CBCT-based adaptive proton planning using an automated planning software for treatment of HN cancer.CBCT-based sCT was demonstrated to be a powerful tool for accurate proton dose calculation in adaptive proton therapy.Contour correction is recommended to reduce the uncertainty for dosimetric parameter estimation, especially for OARs with relatively small volumes.It was found that adaptive IMPT using automated planning of HN cancer resulted in better target coverage and sparing for larynx and constrictor OARs compared to nonadaptation.The NTCP of larynx edema was significantly reduced compared with non-adaptive IMPT, but the magnitude of potential benefits from APT were patient dependent.

Fig. 1
Fig. 1 Workflow for the validation of CBCT-based proton dose calculation on the sCT

Fig. 2 Fig. 3
Fig. 2 An example case with different CT images.a pCT; b rCT; c CBCT; d rCT def with gold standard contours; e sCT with DIR-propagated contours

Fig. 4 Fig. 5 Fig. 6
Fig.4 Box plots of the difference in dose-volume indices between rCT def and pCT, sCT with physician-corrected contours, and sCT with uncorrected contours.The blue box represents the difference between rCT def and pCT, the red box represents the difference between rCT def and sCT with uncorrected contours, and the green box represents the difference between rCT def and sCT with physician-corrected contours

Fig. 7 Fig. 8
Fig. 7 Box plots of NTCP difference between the adapt and non-adapt cohorts

Table 2
Dose constraints for the OARs for IMPT planning

Table 3
Average difference of dose-volume indices between dose calculated on rCT def , pCT, and sCT with physician-corrected and uncorrected contours (range is reported in brackets)

Table 4
Average dose-volume indices and NTCP for planned, non-adapted dose accumulation (non-adapt), and adapted dose accumulation (adapt) A P value < 0.05 is highlighted in bold and indicates statistical significance.A P value < 0.1 is highlighted in italic